JUL  16 


To  be  presented  at  the  Third  Workshop  on  Structural  Health  Monitoring,  Stanford,  CA, 

September  12-14,  2001 


ACTIVE  WAVE  PROPAGATION  AND  SENSING  IN  PLATES 

Anindya  Ghoshal*,  William  N.  Martin\  Maimur  J.  Simdaresan*,  Mark  J.  Schulz^, 

and  Frederick  Ferguson* 


20011022  052 


*NASA  Center  For  Aerospace  Research,  Department  of  Mechanical  Engineering, 
North  Carolina  A&T  State  University,  Greensboro,  NC  27411,  ^Department  of 
Mechanical  Engineering,  University  of  Cincinnati,  Cincinnati,  OH  45221 


D1STR5BUTSON  STATEMEi^T  A 
Approved  for  Public  Release 
Distribution  Unlimited 


ABSTRACT 


Health  monitoring  of  aerospace  structures  can  be  done  using  an  active  interrogation 
approach  with  diagnostic  Lamb  waves.  Piezoelectric  patches  are  often  used  to 
generate  the  waves,  and  it  is  helpful  to  understand  how  these  waves  propagate 
through  a  structure.  To  give  a  basic  understanding  of  the  actual  physical  process  of 
wave  propagation,  a  model  is  developed  to  simulate  asymmetric  wave  propagation 
in  a  panel  and  to  produce  a  movie  of  the  wave  motion.  The  waves  can  be  generated 
using  piezoceramic  patches  of  any  size  or  shape.  The  propagation,  reflection,  and 
interference  of  the  waves  are  represented  in  the  model.  Measuring  the  wave 
propagation  is  the  second  important  aspect  of  damage  detection.  Continuous  sensors 
are  useful  for  measuring  waves  because  of  the  distributed  nature  of  the  sensor  and 
the  wave.  Two  sensor  designs  are  modeled,  and  their  effectiveness  in  measuring 
acoustic  waves  is  studied.  The  simulation  model  developed  is  useful  to  understand 
wave  propagation  and  to  optimize  the  type  of  sensors  that  might  be  used  for  health 
monitoring  of  plate-like  structures. 


INTRODUCTION 

Active  propagation  of  Lamb  waves  in  the  plane  of  the  material  is  a  common 
technique  used  to  detect  damage  [1-7]  in  aerospace  structures.  Methods  of 
measuring  Acoustic  Emissions  (AE)  are  also  used  for  damage  detection  [9-11].  The 
modeling  of  wave  propagation  is  often  done  for  plates  that  are  of  infinite  dimension 
[13-18]  because  closed-form  solutions  of  wave  propagation  in  a  bounded  medium 
are  difficult  to  obtain.  A  normal  mode  approach  for  obtaining  a  closed  form  solution 
for  bending  wave  propagation  in  a  plate  is  presented  in  [19,  20].  Finite-element 
methods  are  also  used  to  model  wave  propagation,  but  the  computation  time  is 
usually  large  and  the  representation  of  higher  frequencies  and  modes  are 
constrained  by  the  number  of  elements  used. 

The  normal  mode  approach  is  also  used  in  the  present  paper  to  model 


asymmetric  Lamb  wave  or  bending  wave  propagation  in  a  simply  supported  plate. 
A  composite  material  plate  is  studied  here,  and  piezoceramic  actuators  and  sensors 
are  modeled  on  the  plate  for  generating  and  receiving  the  waves.  The  objective  of 
developing  this  model  is  to  provide  a  tool  to  simulate  wave  propagation  and  to  help 
the  design  of  sensors  to  measure  the  waves.  The  model  uses  a  closed-form  solution 
and  it  is  written  in  matrix  and  vector  forms  to  run  efficiently  on  a  PC  computer. 
With  this  technique,  series  and  array  configurations  of  sensors  can  be  investigated 
including  different  shapes  of  sensors  and  actuators.  Piezoelectric  patches  are  used  as 
the  sensor  and  actuator  material.  The  sensor  design  can  mimic  receptors  that  excite 
dendrites  that  are  the  inputs  of  neurons  in  the  human  nervous  system.  Ten  or  more 
nodes  can  be  connected  in  a  continuous  single  channel  nerve.  This  continuous 
sensor  approach  is  investigated  first  in  the  paper.  Then,  the  same  sensor  elements 
used  for  the  continuous  sensor  are  connected  in  an  N-by-N  array  that  causes  the 
individual  signals  fi'om  the  sensors  to  be  combined  into  2N-1  array  outputs,  as 
compared  to  the  single  output  from  before.  For  large  array  sizes,  this  approach 
greatly  reduces  the  number  of  channels  of  data  acquisition  instrumentation  needed 
for  structural  health  monitoring.  A  trade-off  in  these  two  approaches  is  that  the 
continuous  sensor  is  the  simplest  with  only  one  channel  of  data  acquisition,  while 
the  array  uses  more  channels  to  more  accurately  locate  damage.  The  modeling  of 
the  plate  and  the  electrical  modeling  of  the  sensors  are  discussed  next. 


MODELING  THE  PLATE 

An  algorithm  to  simulate  waves  propagating  in  a  plate  has  been  developed.  The 
algorithm  uses  a  step  input  to  a  simply  supported  plate,  and  it  takes  into  account 
(integrates)  the  exact  strains  over  the  area  of  the  sensor  node  to  compute  the  voltage 
output.  Details  of  the  plate  and  the  surface  bonded  actuator  model  will  be  given  in  a 
following  journal  paper.  From  the  classical  theory  of  vibration  of  plates,  the 
equation  of  motion  has  been  derived  and  written  in  terms  of  internal  plate  flexural 
moments  Mx,  Mxy,  and  My,  and  actuator  induced  moments  mx  and  my.  The  equation 
of  motion  in  then  expressed  in  displacement  form  as: 
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The  actuator  induced  moments  are  expressed  as  Heaviside  step  functions.  Using 
orthogonalization  and  separation  of  spatial  and  time  variables,  the  temporal 
equation  in  terms  of  modes  is  derived.  The  total  solution  for  the  plate  displacement 
is  expressed  in  terms  of  a  Fourier  series  (Navier  solution)  and  by  summation  of 
modes  for  an  impulse  moment  actuation: 

=  Osin^^)smff)  (2) 

,  .  ffl..  "  «  b 


where  o),  =  o)„„  yjl-CL  ,  and 


F  mn  — 


-r-r  (.-^)['»x  —  +  my  — ](cos( - '- 

phab  na  mb  a 


)  -  cos(-^^))(cos(-^)  -  cos(-^^)) 
abb 


(3) 


where  a  and  b  are  the  width  and  breadth  of  the  plate  in  x  and  y  directions,  m  and  n 
are  the  mode  numbers,  h  is  the  thickness  of  the  plate,  m\  and  m\,  are  the  distributed 
surface  moments,  p  is  the  mass  density  of  the  fiberglass  plate,  emd  z  is  h/2.  The 
corresponding  strains  are: 
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MODELING  OF  THE  SENSOR 

The  piezoceramic  sensor  nodes  can  be  modeled  as  a  capacitor  in  parallel  with  a 
current  source.  The  piezoelectric  constitutive  equations  used  are  listed  in  the  IEEE 
standard  ANSI/IEEE  Std.  176-1987.  These  equations  can  be  put  into  matrix  form  to 
give: 
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The  piezoceramic  sensors  can  be  modeled  using  the  piezoelectric  constitutive 
equations  and  by  connecting  the  segments  into  an  electric  circuit  as  shown  in  Fig.  1. 
Please  refer  to  [12]  for  the  full  derivation,  but  the  expression  for  the  output  voltage 
equation  for  a  sensor  nerve  is: 
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where  C  is  the  capacitance  of  the  PZT  sensor,  the  effective  capacitor  area  is  Ag,  the 
effective  plate  separation  distance  is  h,  ic  represents  the  component  of  the  current 
going  through  the  capacitor  of  the  model,  and  ig  represents  the  component  of  the 
current  generated  by  the  piezoelectric  sensor.  The  homogeneous  and  particular 
solutions  for  Equation  (9)  must  be  calculated  and  added  for  the  total  solution  of  the 
current  i.  The  product  of  the  current  i(0,  and  the  impedance  R  of  the  measuring 
device  equals  the  voltage  of  the  series  connected  sensors  as  a  function  of  time.  Thus 


we  solve  for  the  current  to  get  the  voltage  V^  =  \R.  This  voltage  is  proportional  to  the 

dynamic  strain  in  the  structure  at  the  sensor  and  thus  can  be  used  to  detect  damage 
through  dynamic  strain  measurements  and  acoustic  emissions. 


Figure  1.  Circuit  model  of  sensor  (a)  continuous  sensor,  (b)  cross-sensor  array. 


SIMULATION  OF  WAVE  PROPAGATION 

The  simulation  is  performed  using  a  model  of  the  fiberglass  plate.  Three  different 
cases  of  active  wave  propagation  using  a  surface  bonded  PZT  actuator  and  sensing 
using  continuous  and  array  sensors  are  presented  here.  The  material  properties  of 
the  fiberglass  plate  can  be  found  in  [1 1].  The  size  of  the  plate  is  0.88m  x  1.21m  x 
3.2mm.  The  first  100  vibration  modes  have  been  used  for  the  simulation  and  the 
time  step  used  is  one  micro-second. 

Figure  2  shows  the  active  wave  propagation  at  the  times  10  micro-sec  to  360 
micro-sec  due  to  anti-symmetric  Lamb  waves  generated  by  a  PZT  actuator  placed  at 
the  center  of  the  fiberglass  plate.  The  actuator  dimensions  are  2.5cm  x  5cm  x 
0.25nim.  Figure  3  shows  the  voltage-time  history  of  the  continuous  and  the  array 
sensors  modeled  on  the  plate. 

Figure  4  shows  the  active  wave  propagation  at  the  times  10  micro-sec  to  360 
micro-sec,  due  to  anti-symmetric  Lamb  waves  generated  the  PZT  actuator  placed  at 
one  edge  of  the  fiberglass  plate.  The  actuator  dimensions  are  2.5cm  x  5cm  x 
0.25mm.  Figure  5  shows  the  voltage-time  history  of  the  continuous  and  the  array 
sensors  modeled  on  the  plate. 

Figure  6  shows  the  active  wave  propagation  at  the  times  10  micro-sec  to  360 
micro-sec,  due  to  anti-symmetric  Lamb  waves  generated  the  PZT  actuator  placed  at 
the  center  and  across  the  width  of  the  fiberglass  plate.  The  actuator  dimensions  are 
0.81m  X  2.5cm  x  0.25mm.  Figure  7  shows  the  voltage-time  history  of  the 
continuous  and  the  array  sensors  modeled  on  the  plate. 

Currently  efforts  are  imderway  to  model  a  crack  in  the  plate  and  use  the  actuator 
and  sensor  bonded  on  the  plate  to  detect  the  crack.  Further  investigations  are  being 
done  to  correlate  the  peaks  in  the  voltage  time  history  graphs  with  the  time  when 
the  incident  and  the  reflected  waves  hit  the  sensors,  and  when  the  waves  are 
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scattered  from  a  crack  embedded  in  the  plate.  A  neural  network  algorithm  is 
currently  being  devised  to  detect  and  locate  the  erack  in  the  plate  model. 


Figure  2.  Wave  propagation  at  10, 60, 110, 160, 210, 260, 310,  360  micro-sec  due 

to  center  actuation  by  a  PZT  patch. 


Figure  3.  Voltage  time  history  due  to  center  aetuation  by  a  PZT  patch;  (a)  array 
voltage  for  2  ms,  (b)  array  voltage  for  0.2  ms. 


Fig  5.Wave  propagation  at  10,  60,  1 10, 160,  210, 260,  310,  360  micro-sec  due  to 

edge  actuation  by  a  PZT  patch. 
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Figure  6.  Voltage  time  history  due  to  edge  actuati 
voltage  for  2  ms,  (b)  array  voltage 
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igure  7.  Voltage  time  history  due  to  edge  actuation 
sensor  voltage  for  2  ms,  (b)  continuous  sens 


Figure  8,  frames  1-4 


Figure  8  continued.  Wave  propagation  at  10,  60,  110, 160, 210, 260,  310, 360 
micro-sec  due  to  center  actuation  by  a  wide  PZT  patch. 
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Figure  9.  Voltage  time  history  due  to  center  actuation  by  a  wide  PZT  patch;  (a) 
array  voltage  for  2  ms,  (b)  array  voltage  for  0.2  ms. 


Figure  10.  Voltage  time  history  due  to  center  actuation  by  a  wide  PZT  patch;  (a) 
continuous  sensor  voltage  for  2  ms,  (b)  continuous  sensor  voltage  for  0.2  ms. 


CONCLUSIONS 


The  simulations  perfonned  have  shown  that  multiple  piezoceramic  patches  can  be 
connected  together  in  a  series  or  array  pattern  to  simulate  the  way  biological  nerves 
have  multiple  inputs  (dendrites)  cormected  together.  This  reduces  the  number  of 
channels  of  data  acquisition  needed  to  detect  damage  represented  by  acoustic 
emissions  or  high  strains.  The  simulation  model  developed  can  be  used  to  optimize 
the  design  of  the  neural  system  for  different  structural  materials  and  sizes.  This  can 
be  used  for  simulating  active  interrogation  for  damage  detection  in  plate  type 
structures.  It  helps  in  optimizing  the  actuator  and  sensor  locations  for  both  the 
continuous  sensor  nodes  and  for  the  sensor  array  configurations.  The  model  can 
also  be  used  for  devising  a  neural  network  algorithm  for  detection,  localization,  and 
quantification  of  the  crack/damage  and  monitoring  the  survivability  of  the  plate 
structure. 
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